Entanglement and entangling power of the dynamics in light-harvesting complexes 
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We study the evolution of quantum entanglement during exciton energy transfer (EET) in a network model of 
the Fenna-Matthews-Olson (FMO) complex, a biological pigment-protein complex involved in the early steps 
of photosynthesis in sulphur bacteria. The influence of Markovian, as well as spatially and temporally correlated 
(non-Markovian) noise on the generation of entanglement across distinct chromophores {site entanglement) and 
different excitonic eigenstates (mode entanglement) is studied for different injection mechanisms, including 
thermal and coherent laser excitation. Additionally, we study the entangling power of the FMO complex under 
natural operating conditions. While quantum information processing tends to favor maximal entanglement, near 
unit EET is achieved as the result of an intricate interplay between coherent and noisy processes where the initial 
part of the evolution displays intermediate values of both forms of entanglement. 



I. INTRODUCTION 

Photosynthesis, at its simplest, is the absorbtion of sunlight 
by photosensitive antennae and its subsequent conversion into 
chemical energy at a reaction center. The locations for these 
processes are physically and physiologically separated, which 
forces nature to devise a way for transferring the solar energy 
from the antennae to the reaction center (RC). Exciton en- 
ergy transfer (EET) is facilitated by certain protein molecules 
called light-harvesting complexes, and occurs at efficiencies 
of about 99%. EET has been a subject of continual interest for 
decades, not only for its phenomenal efficiency but also for 
its fundamental role in Nature iQl] . Recently, ultrafast optics 
and nonlinear spectroscopy experiments have provided new 
insights into the process of EET in light-harvesting complexes 
like the one found in purple bacteria (LH-I) and the Fenna- 
Matthew-Olson (FMO) complex In particular, evidence 

of quantum coherence has been presented, with the idea that 
nontrivial quantum effects may be at the root of its remarkable 
efficiency |4[|. Following this, several studies have attempted 
to unravel the precise role of quantum coherence in the EET of 
light-harvesting complexes Ipi- Ulll . and have, perhaps surpris- 
ingly, found that environmental decoherence and noise plays 
a crucial role BMflUl. 

Light-harvesting complexes consist of several chro- 
mophores mutually coupled by dipolar interactions residing 
within a protein scaffold. Due to their mutual coupling, 
light-induced excitations on individual chromophores (sites) 
can undergo coherent transfer from site to site, and the typ- 
ical eigenstates are therefore delocalized over multiple chro- 
mophores. It is in these eigenstates, henceforth referred to 
as exciton states, that one finds evidence of quantum coher- 
ence. Here, we will study the role of quantum coherence in 
the process of EET in the FMO complex, as quantified by 
quantum entanglement, and investigate the sensitivity of the 
entanglement dynamics to variables that have not been di- 
rectly measured, such as the microscopic interaction strength 
between the complex and the surrounding environment, and 
the possible existence of spatial and temporal correlations in 
the bath. In this context, the first analysis of the entangle- 



ment behaviour in light-harvesting complexes was presented 
in Ref. in which we analyzed the evolution of an entan- 
glement measurement, i.e. logarithmic negativity lfl3i[l4ll . in 
a Markovian model of the FMO complex. The scope of the 
present work is two-fold, on the one hand it aims to make the 
study of coherence and entanglement in such systems quan- 
titative by considering spatial and temporal (non-Markovian) 
noise correlations, and secondly, it uses this quantitative ap- 
proach to show that maximal entanglement is not correlated 
with optimal transport, a result that may shed light on the pos- 
sible functional role of entanglement in EET. 

Entanglement is defined between subsystems of a global 
system. When considering entanglement in a composite sys- 
tem whose components are closely spaced and strongly inter- 
acting, as in the FMO complex, this choice of subsystems is, 
to some extent, dictated by the way we interrogate the sys- 
tem. If the sites can be addressed individually, then it is oper- 
ationally well justified to speak about site-entanglement, i.e. 
quantum correlations across distinguishable locations. How- 
ever, if we are limited to accessing the global excitations of 
the systems, i.e. the excitonic eigenstates of the Hamilto- 
nian governing the dynamics of the FMO complex, then we 
will speak of mode entanglement as we then explore entan- 
glement between the eigenmodes of the system (for more de- 
tails on entanglement theory and its essential foundations see 
Ref. ifTill ). Here, we will explore both types of correlations 
but place a perhaps greater emphasis on the aspects of site- 
entanglement which is more closely linked to the non-local 
structure of quantum correlations. In general, the presence of 
quantum coherence, signified by the presence of off-diagonal 
elements in the density matrix, is necessary but not sufficient 
for the presence of quantum entanglement 13 ■ However, the 
two conditions are equivalent if one makes the idealized as- 
sumption that there is a single excitation in the system. In this 
work, as explained below, our entanglement analysis goes also 
beyond the single-excitation approximation and takes into ac- 
count different experimental and natural operating conditions 
for the FMO complex dynamics. 

The paper is organized as follows. In Sec. |ll] we intro- 
duce the theoretical model for the FMO dynamics and discuss 
the entanglement measure used here. Then, we analyze the 
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logarithmic-negativity and the so-called entangling power in 
some non-Markovian FMO models in Sec. |III]and generalize 
these results for different injection schemes, i.e. thermal in- 
jection and laser excitation of the FMO complex, in Sec. HV] 
Finally, the entanglement between excitons (mode entangle- 
ment) is investigated in Sec. |V]and the conclusions and final 
remarks are presented in Sec. |Vll 

11. SITE ENTANGLEMENT 

Each chromophore (site) in the FMO complex can be rep- 
resented by a two-state system (qubit). As in Ref. 01, the 
effective dynamics is modelled by an = 7 qubit Hamil- 
tonian which describes the coherent exchange of excitations 
between sites, i.e., 

7 

H = fiwjcr+cr^ + hVj^i{(J-(J+ + cr+crf ) , (1) 

and local Lindblad terms that take into account the dissipation 
and dephasing caused by the surrounding environment, i.e., 

7 

Cdtsa{p) =Y^][-{(^t^J + '^^JP^t^ (2) 

and 

7 

Cdephip) ^^iji-i'^t'^J 'P^ + '^'^t'^7 P'^t'^7^ ■ 

This Markovian form of the evolution preserves complete pos- 
itivity, an essential feature when evaluating entanglement, as 
discussed below. Here cr^ i'^J) the raising and lowering 
operators for site j, hutj is the local site excitation energy, Vk.i 
denotes the hopping rate of an excitation between the sites k 
and I, and Tj and 7^ are the dissipative and dephasing rate at 
the site j, respectively. Finally, the transfer efficiency is mea- 
sured in terms of an irreversible transfer of excitations (with 
rate Tsink — 6.283 ps~^) from site 3 to an extra site 8 model- 
ing the RC, as described by the Lindblad operator 

Particularly, the transport efficiency is described by the popu- 
lation transferred to the sirLkpsmfc(i), which is given by 

Psink{t) ^2Tsink I Pz{t')At' , (5) 

where p-i{t') is the population of the site 3 at time t' . More- 
over, we choose Tj = 5 x 10"'* ps~^ for any site j, as in 
Ref. Notice that decoherence appears in the model above 
via the action of a pure-dephasing Lindblad super-operator in 
the master equation for the exciton dynamics, which is equiv- 
alent to having stochastic fluctuations of the exciton site en- 
ergies induced by the environment. This model is known as 
the Haken-Strobl model and has been used extensively in the 



chemical physics literature to describe exciton dynamics over 
several decades 11511 . In the next section, we will investigate 
the entanglement behaviour in some non-Markovian models 
of the FMO complex dynamics with different types of envi- 
ronmental interactions. 

Let us stress that, in order to obtain sensible and reliable 
results for the entanglement analysis, it is crucial that the ex- 
citon dynamics is represented by a completely positive map. 
Indeed, we numerically found that small deviations from the 
complete positivity conditions are enough to cause significant 
changes to the entanglement. The Lindblad formalism abso- 
lutely guarantees the positivity of the evolving state, whereas 
the majority of non-perturbative and non-Markovian treat- 
ments do not. Hence, in the following, we deem it prudent to 
choose and to analyze only non-Markovian models, that give 
both reliable entanglement predictions and which are consis- 
tent with the essential experimental transport data. Actually, 
very little is actually known about the microscopic details of 
the environment in FMO complex, and it is very hard to distin- 
guish between the various noise models proposed in literature 
- see for instance Refs. Jsl [13, [3 [i3l - using the available 
experimental data. One may conclude that the additional com- 
plexity of all of these models results from the introduction of 
a range of new variables which can be all independently tuned 
to match the key experimental results, even though the de- 
tailed dynamics of these models may be quite different. It 
should be noted in this respect that the simple Lindbladian 
model we use here can in fact account for the key features of 
the experimentally observed dynamics, i.e. long-lived coher- 
ence times and transport times - see Ref. 

We will quantify the entanglement across a bipartition A\B 
of a composite system by using the logarithmic negativity 
llll,i.e. 

=log2||/-''||i , (6) 

where Ta is the partial transpose operation of the density op- 
erator p with respect to the subsystem A and 1 1 • 1 1 1 denotes 
the trace norm. It quantifies how negative the spectrum of the 
partial transpose of the density matrix is, consequently it is 
only meaningful if the evolution is completely positive. Fur- 
thermore, we would like to stress that in addition to its com- 
putational simplicity, the logarithmic negativity also possesses 
an operational interpretation in terms of the entanglement cost 
for the exact preparation of the state ifisll . When confined to 
the one excitation subspace, \f A ~ 1 ... fc is a set of fc chro- 
mophores within a global system of N sites, then the logarith- 
mic negativity across the bipartition (1 . . . fc) | (fc + 1 . . . N) is 
given by the compact expression 

. . . fc|fc + 1 . . . iV) = log2 - aoo + ^al^+^X^ 

(7) 

where X ~ Sj=fc+i ky Pi ™d atj denotes the off- 

diagonal element between states with excitations in qubits i 
and j. Here aoo is the matrix element corresponding to the 
zero excitation subspace. If all coherences are vanishing, i.e. 

= 0, there is no entanglement across any partition in the 
one-excitation sector Note that the restriction to at most a sin- 
gle excitation is unproblematic in the case of the logarithmic 
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negativity as it does not affect its definition and the fact that it 
is an entanglement monotone, that is non-increasing, for gen- 
eral local operations and classical communication (LOCC). 
This is not the case if the constraint to at most a single ex- 
citation is applied in a way that amends the definition of the 
functional. In Ref. fl9ll . for example, the global entanglement 
Eg, which is defined as the relative entropy of entanglement 
with respect to the set S of totally separable states 1I20I1 . was 
amended to yield Ed by replacing the set S of totally separa- 
ble states by the set Si of totally separable state with at most 
one excitation. This new function Eqi is non-increasing only 
when we restrict attention to the set of LOCC operations Oi 
that map Si into itself. The set Oi excludes a wide range of 
important physical processes such as non-diagonal local uni- 
taries and hence fundamental processes such as laser excita- 
tion because the raising operator ^ Oi does not map Si 
onto itself. Hence, while being computable Eqi is not an en- 
tanglement monotone under natural, physically realizable and 
important operations that are routinely applied. In fact, the op- 
erations permitted in Oi do not create coherence between the 
zero and the single excitation sector and hence a more natural 
interpretation of Eqi is one of quantifying coherence rather 
than entanglement. 

Therefore we will restrict attention to the logarithmic neg- 
ativity, which is an accepted and at the same computable en- 
tanglement monotone lll4ll for arbitrary LOCC and arbitrary 
excitation levels. 



III. NON-MARKOVIAN MODELS 

To estimate the impact of non-Markovian effects on the en- 
tanglement in the FMO complex, we will consider two non- 
Markovian models of the FMO complex dynamics, corre- 
sponding to a different type of interaction with the surround- 
ing environmental. In particular, we assume the FMO com- 
plex to be linearly coupled to a bath of damped harmonic os- 
cillators. In all cases, the harmonic oscillators will be damped 
into a zero temperature bath at rate k. By using these non- 
Markovian models, which manifestly preserve the completely 
positivity of the corresponding quantum evolution, we can 
study the behaviour of the entanglement, measured by the log- 
arithmic negativity and by the entangling power of the quan- 
tum evolution itself. 



A. Local bath model 

While the phenomenon of BET is clearly noise-assisted, the 
exact dynamics within the observed exciton transmission time 
are strongly model-dependent. To further emphasize this as- 
pect and therefore invoking the need of further experimental 
results, we will present in this section a tunable noise model 
which reproduces the Markovian results presented in |7] in 
certain parameter regime but, interestingly enough, can also 
provide longer coherence times while preserving and even 
enhancing EET in a parameter regime where the model ex- 
hibits a degree of non-Markovianity. In particular, we con- 



sider an environment model motivated by the approach pre- 
sented by Adolphs and Renger in Ref. [30, where sites interact 
locally with a quasi-resonant localized mode. Their spectral 
density contains a contribution from a low-energy continu- 
ous density of states and a discrete high-energy mode, and 
its effects on the dynamics of a dimer molecule were recently 
simulated using a new application of the time-adaptive renor- 
malisation group method |21]. Here, for simplicity, we will 
consider a model in which each FMO chromophore is linearly 
coupled to a resonant harmonic mode with strength g while 
each mode is damped into a zero temperature Bosonic reser- 
voir with strength k. In order to describe these couplings, we 
add to the previous Hamiltonian in Eq. ([T]) the following two 
terms 

7 
7 

HsB = a J {a J + a+ ) (7+ G- , (9) 

where Rb is the free Hamiltonian for the two-level bath with 
creation and annihilation operators a+ and a, respectively, and 
mode frequency = ujj, and Hsb is the system-bath inter- 
action Hamiltonian with interaction strength g. The damping 
is introduced by considering a Lindblad term C}^^^^^ (p) of 
the form 

7 

l^\-aT\p) = E P} + 2a,pa+] (10) 

with kj being the rate at which the local harmonic mode, 
coupled to the site j, is damped into a zero temperature bath. 
Within the considered parameter regimes, the local modes can 
be reasonably considered within a two-level approximation. 
In fact, we numerically monitor the populations in each local 
mode and we found that no local modes were strongly excited 
or saturated over the whole time interval investigated here. 
The damping rate Ki determines the width of the spectrum and 
hence the correlation time of the local environment associated 
to the site %. For weak coupling and strong mode-losses, this 
approach leads to a Markovian environment with Lorentzian 
line shape (similar to the models studied by different methods 
in ll2ll I22I1 ). while, for low losses and strong coupling, 
the high degree of excitation of the environment leads to 
deviations from the Lorentzian lineshape. To isolate the 
impact of the non-Markovianity, we keep the ratio j k., i.e. 
the effective coupling strength between the site and its mode, 
fixed while varying the ratio gj k. To this end, we employ 
a parameter / to parametrize the system-mode coupling 
rates g = \/7.9o and energy loss rates of the modes with 
K = /ko. Then for / 3> 1 (5 ^ k), i.e. the Markovian limit, 
we reproduce the optimized dephasing rates found in |01, 
while for / <C 1 (,9 ^ k) we find non-Markovian behaviour. 
We initiate the system with one excitation in site \. For a 
single site this dynamics has been tested to reproduce both 
the correct Markovian limit and to be capable of exhibiting 
strongly non-Markovian behaviour as illustrated in Fig. 
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FIG. 1. Entanglement of different bi-partitions of the FMO com- 
plex dynamics subject to a local form of noise where each FMO 
site couples to a damped resonant mode. The effective dephasing 
rates are taken to be equal to the optimal values derived for a fully 
Markovian model as in Ref. |0] but the degree of Markovianity of 
the model can be tuned by varying the parameter /, as described in 
the text. The entanglement curves are for the 6 splits of the form 
(1, • • • , + 1, • • • , 7), i £ {1, • • ■ , 6}, where i corresponds to 
the FMO site i. The behaviour of Psi„k for different values of / is 
also shown. Large values of / render the dynamics fully Markovian 
and we recover the results in Ref. 0], while one can observe longer 
coherence times while at the same time enhancing Psink in the low 
/ domain, as the one exemplified by the value / = 1. Note the 
non-monotonic behaviour of the transfer rate as a function of /, as 
emphasized by the arrows. 
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FIG. 2. Entanglement of different bi-partitions of the FMO complex 
vs. time (ps), in the presence of a tunable source of noise correla- 
tions described by a non-local Bosonic bath, for different values of 
the ratio /. For small /, the environmental correlations yield a non- 
Markovian form of noise and the time scale for entanglement persis- 
tence increases while, within this model, the efficiency for transport 
Psink decreases. 



of the microscopic details of the environment. In this respect, 
experimental entanglement measurement could be important 
for removing some of this uncertainty, as it can be very 
sensitive to the structure of the noise the system is being 
subject to ll24ll . 



[U Here, we choose the system-mode coupling rates to be 
equal to 50 = = {1, 50, 41, 50, 41, 5, 50}/5.3 ps^^ 
to match closely the effective dephasing rates 
{0.157,9.432,7.797,9.432,7.797,0.922,9.433} ps^^ in 
the Markovian limit, as in Ref. H. In Fig. □ we show 
the entanglement behaviour for two different bi-partitions 
of the FMO complex and for two values of / far from strict 
Markovianity as well as the time-dependence for Psink when 
/ varies from / = 0.1 to the value of 100. The behaviour 
obtained in the fully Markovian case of Ref. ^ is also 
represented and is indeed already reached for / = 100. Inter- 
estingly, we obtain that taking / as our measure of deviation 
from Markovian behaviour, Psink varies non-monotonically 
with / so that, differently from the non-local bath non- 
Markovian model below or the model described in ll22ll . here 
the presence of a degree of temporal correlations may even 
assist the transport of electronic excitations from the antenna 
to the RC, as seen for values of / in the range 1 — 10. On the 
other hand, very low values of / (strong non-Markovianity) 
lead to decreased transport while preserving large values of 
the coherence across site bipartitions. Note that different 
noise models, as the one presented in ll23ll . yield different 
conclusions. This also shows the great uncertainty about the 
nature of transport dynamics that arises from our ignorance 



B. Non-local bath model 

In this subsection we would like to go further and explore 
more sophisticated entanglement properties of the entangle- 
ment dynamics including its entangling power To this end 
we need to further simplify the decoherence model to permit 
its numerical analysis in the context of the entangling power. 
This model will also include non-local correlations in the en- 
vironment and we will start by considering again the entan- 
glement of states in the FMO dynamics under this model and 
then apply it to the study of the concept of entangling power. 

Formally, we consider a model in which we couple each site 
linearly, with strength gj, to a single, common harmonic mode 
damped into a zero temperature bath at rate Kj that depends on 
the location of the excitation ll25ll . In particular, the following 
bath-system interaction Hamiltonian term Hsb is added to 
previous Hamiltonian in Eq. ([U, i.e. 

7 

Hsb = 9o [a + a+)a+aj , (11) 
and the damping is described by the following Lindbladian 
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time (ps) 

FIG. 3. Contour plot for the transfer efficiency in FMO complex, as a 
function of time (in ps) and of the parameter /, in the case of a non- 
Markovian model with a non-local Bosonic bath linearly interacting 
with the FMO. 



super-operator 

7 

-C^a'^^P) = Y.K,[-{¥,a+a¥,,p} + 2¥,apa+¥,] , (12) 

where a and a+ are, respectively, the annihilation and cre- 
ation operators of the harmonic mode and Vj is the projec- 
tor onto site j in the FMO complex. Hence the line width 
of the harmonic oscillator will depend on the location of the 
excitation. This model mimics closely the model of the pre- 
vious section in that sites see harmonic oscillators with site- 
dependent damping rates. At the same time, the Hilbert space 
of this model is much smaller as the dynamics is restricted 
to the single excitation sector only. In effect we can describe 
the dynamics in the basis {|i)|0), . . . , |i) |(i)}i=i...._7 where the 
first index refers to the excitation in the FMO complex and 
the second refers to the environment oscillators. The damping 
rate of the environment oscillator can be made to depend on 
the site in which the excitation resides in as would be the case 
in the local model. Hence this model represents a mix be- 
tween the local model in the previous subsection and that of 
an FMO complex coupled to a single mode the former allow- 
ing for different coupling rates and line widths for different 
sites and the latter allowing for a significant reduction in the 
Hilbert space dimension of the simulation. In Fig. [2] we con- 
sider the entanglement between site 1 and the remainder of the 
FMO complex as well as sites 1 and 2 versus the remainder of 
the FMO complex. Here, we choose the system-mode cou- 
pling rates as in the local non-Markovian model above, and 
we find that non-Markovian effects (decreasing /) reduces the 
transport efficiency while it prolongs the lifetime of entangle- 
ment. These observations are explained by the fact that the 
non-Markovian dephasing leads to a reduction of the effec- 
tive noise level in the system (a phase flip by the environment 
may be followed by another correlated phase flip at a later 
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FIG. 4. Contour plot for the entangling power in FMO complex, 
as a function of time (in ps) and of the parameter /, in the case 
of the non-local bath non-Markovian model. For low noise levels 
the entangling power of the evolution can be non-monotonic in time, 
indicating coherent oscillations in the system. The entangling power 
decreases monotonically with increasing noise level. As comparison, 
the entangling power of a CNOT gate is 1, and a SWAP gate is 2. This 
means that the FMO complex has nontrivial entangling power for a 
large parameter range. 



time, hence canceling out), upsetting the optimal balance of 
quantum and incoherent dynamics required for efficient EET 
Ei^, whilst also increasing and preserving the entanglement 
that is present in the system. The fact that the presence of a 
non-Markovian environment may enhance entanglement con- 
tent beyond the values predicted for an evolution subject to 
memoryless environments was also predicted when analyzing 
strictly bipartite systems ll26l l27ll . including biological sce- 
narios ll28ll . Moreover, rather than exploring the entanglement 
content of states that are generated during the evolution, one 
may also study quantitatively the entanglement content of the 
quantum evolution itself. Just as for quantum states such a 
quantification is, however, not unique. Here, we consider the 
capacity of an evolution to create entanglement between two 
subsystems, each of which may be composed of several com- 
ponents. This automatically provides a lower bound for the 
amount of entanglement that is required to reproduce the dy- 
namics of the syster n ptirely from local operations and classi- 
cal communication ll29ll3C)ll . In Figs.[3]and|4] we show a con- 
tour plot for the transfer efficiency and the entangling power 
of the noisy evolution of the FMO complex in the first picosec- 
onds of the EET process, respectively. In particular, we con- 
sider the entangling power for the evolution between different 
bi-partitions of the system, quantified by the logarithmic neg- 
ativity, in the FMO when initially prepared in a maximally 
entangled state with 7 ancilla qubits. We consider the split 
{qubits lFA/o-la?icjna}-{the rest}, as a function of time (in 
ps) and of the parameter above /. Large quantum correlations 
are not associated with optimal transport. In the absence of 
any dephasing, entanglement lasts for times limited only by 
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the excitation loss rate. On the other hand, non-Markovian 
dephasing (decreasing /) increases the entanglement, which 
persists for about the initial 20% of the total transmission time, 
but decreases transport efficiency. 



IV. BEYOND SINGLE EXCITONS 

So far we have assumed that the system is initialised with a 
single excitation in site 1 . This may not be realised precisely 
under experimental or natural operating conditions and, fur- 
thermore, neglecting higher excitations (though only existing 
for a really short time) may influence the entanglement con- 
tent of the system considerably. Hence, a study of quantum 
entanglement in the FMO complex under realistic conditions 
should consider a model which allows the freedom to control 
the number of excitations in the complex at any time. To this 
end, we generalize the theoretical noise model proposed in 
ialTtl, by first modeling (i) the baseplate feeding excitations 
into the FMO complex as a thermal reservoir of excitations at 
an effective temperature T and then (ii) a system under laser 
pulse irradiation. 




time (ps) 

FIG. 5. Entanglement, quantified by the logarithmic negativity, in 
the FMO complex for local (dashed lines) and for spatially corre- 
lated (continuous lines) dephasing noise, in presence of the thermal 
injection bath (nth ~ 100, Fi = 1), within the Markov approxima- 
tion. The inset shows the site population behaviour as a function of 
time. 



A. FMO thermal injection 

The complex starts in the ground state, without any exci- 
tations, which are introduced into the network via the site i 
with a rate F^. This process is modelled by a thermal bath of 
harmonic oscillators at a temperature given by the thermal av- 
erage boson number nth- Within the Markov approximation, 
the Lindblad superoperator for the injection of excitations as- 
sumes the form 

p 

+ (nth + l)y + 2arpa+] , (13) 

which can now be used to study the evolution of quantum 
entanglement in the FMO complex (or any light-harvesting 
complex) under various possible natural settings. In partic- 
ular, motivated by the experimental observations, we choose 
the site 1 as the one at which the excitations are introduced 
in the FMO complex. In Fig. |5] we show the entanglement 
time evolution for various subsystems of the qubit network 
modelling the FMO complex. As can be seen, the amount of 
entanglement is considerably smaller compared to when the 
FMO complex starts with exactly one excitation, say on site 1, 
although it persists on the same timescale This is because 
the second term in the injection Liouvillian, which allows for 
an irreversible loss of excitations from the complex (compa- 
rable to spontaneous decay), always accompanies the first one 
which introduces excitations into the complex. Besides, given 
the dimensions and structure of the FMO complex (closest 
site-site distance is ~' 11 A), dephasing may not be local, but 
correlated in space and time la [231 ■ For this reason, 

we consider also the case of spatially correlated dephasing by 



using the following Lindblad term: 

^dephip) = - 5l7nm[An, K, p]] , (14) 

where 7 is a positive semidefinite matrix (to preserve the com- 
plete positivity of the quantum evolution), but with the diag- 
onal elements equal to the optimal local dephasing rates as in 
Ref. 101, and A,n — cFm^m- Here, the amount and duration of 
entanglement are both slightly enhanced (solid lines in Fig.|5]). 

B. FMO laser excitation 

In the laboratory 10, 01, the complex is typically irradiated 
with a short laser pulse centered on the typical transition fre- 
quencies of the sites. The coupling between the FMO com- 
plex and the external radiation field can be described by the 
semiclassical time-dependent Hamiltonian Hp^fo-iaser{t), 
which reads in rotating wave approximation as follows 

7 

i/FA/0-iaser(t) =-5Z^'^eS(t)e-^"i*a+ + /l.C. (15) 
1=1 

with jli being the molecular transition dipole moment of the 
individual site i (taken from the published crystal structure in 
Ref. 1I33I1 ). e being the polarization of the field, and E{t) the 
time-dependent electric field. In particular, following Ref. IH, 
we consider a Gaussian electric field pulse of width 60 fs, cen- 
tered at 120 fs, with an electric field strength Eq = 4.97968 
cm~^, polarized parallel to the dipole moment of site 
1 , and with a frequency on resonance with the optical transi- 
tion of site 1 (cji). The electric field amplitude was chosen to 
excite one excitation on site 1, i.e. to give a tt -pulse on site 
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FIG. 6. Entanglement in FMO complex (evolving by a Markovian 
dynamics) when irradiated with a laser pulse (see the text for more 
details). Dephasing noise is either local (dashed lines) or spatially 
correlated (continuous lines). Strong entanglement between sites 1 
and 2 is clearly illustrated by the behaviour observed in the top curve 
(bipartition 1 — (2 — 7)) as opposed to the bottom one (bipartition 
(1 — 2) — (3 — 7)). Similarly, strong entanglement between sites 
5 and 6 is best exemplified by the central plot where we compute 
the entanglement across bipartition (1 — 5) — (6 — 7). Inset: site 
population distribution in the FMO complex for t = 75 fs. 



1. However, as the molecular transition dipole moments are 
of the order of the intersite energy difference, the laser pulse 
leads to excitation of all sites. We have studied this scenario 
as well, both for local and correlated spatial dephasing noise, 
and the results are presented in Fig. |6] This scenario gen- 
erates a considerable amount of entanglement (though for a 
shorter amount of time compared to the previous scenario), 
and it might be concluded that although an FMO complex 
operating in nature may not possess substantial amounts of 
entanglement, it is possible in a laboratory to generate large 
amounts of it. This could open up new vistas for exploration 
of quantum effects in biological systems, albeit under labora- 
tory conditions, and also allow for a demonstration of entan- 
glement enhancement under non-Markovianity. 



V. MODE ENTANGLEMENT 

When local addressing is unfeasible, site entanglement is 
directly immeasurable even when the evolution is confined 
to the single excitation sector. However, the evaluation of 
mode entanglement via the experimental determination of co- 
herences in the exciton basis can provide us with informa- 
tion about the existence of quantum correlations in the sys- 
tem. The temporal behaviour of mode entanglement within 
our Markovian model 0] and for different bi-partitions is 
shown in Fig. |7] together with the exciton populations along 
the first ps of the BET Initially all modes are populated, with 
the largest fraction in excitons 3 and 6, which leads to high 
values of mode entanglement across bipartitions each contain- 
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FIG. 7. The logarithmic negativity in the FMO complex (evolving 
by a Markovian dynamics) for local dephasing noise in the exciton 
basis. Initially, one excitation is in the site 1. The curves are for 
the 6 splits of the form (1, ■ • ■ , + 1, • • ■ ,7), i £ {1, ■ • ■ , 6}, 
where now i corresponds to the exciton i and they are ordered with 
increasing exciton energies. Inset: exciton population vs. time (ps). 



ing one of those high energy excitons. As time elapses, en- 
tanglement degrades monotonically as the transfer efficiency 
increases. 



VI. CONCLUSIONS AND OUTLOOK 

Efficient EET in light-harvesting complexes can be traced 
back to an interplay between coherent and incoherent pro- 
cesses where the quantum correlations characteristic of the 
coherent evolution are partially suppressed by noise, yet not 
entirely destroyed. We have placed the analysis of entangle- 
ment in such systems on a quantitative footing and showed 
that for optimal transport the entanglement, while present, is 
neither maximal nor long lived. Actually, long lived entan- 
glement exists in the absence of dephasing which is known 
to be highly inefficient. However, despite the fact that ob- 
served transfer times do require a noise-assisted transport dy- 
namics, it turns out that the time-dependence of both coher- 
ences and population transfer over the full transfer time are 
strongly model-dependent. 

In summary, an interplay between creation of entanglement 
for short distances and times (through coherent interaction) 
followed by the destruction of entanglement for longer dis- 
tance and times (through dephasing noise) seems to be nec- 
essary for optimal transport. Moreover, our entanglement re- 
sults could actually be seen as providing a potential experi- 
mental test for the form of the system-environment coupling. 
Unlike the coherence and transport times, the entanglement 
is sensitive to the precise evolution of the system, and an ex- 
periment that could measure entanglement (such experiments 
are indeed planned), could differentiate between models. In 
order to do so, the entanglement predictions for various noise 
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models must be available for comparison, and in this respect 
the nature of the environment is central to the problem of en- 
ergy transport in photosynthetic complexes. Further studies 
are however required before any result can be accepted as 
conclusive on the functional and possibly beneficial role of 
coherence and entanglement in EET. 
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